Quantum Electrodynamics of Internal Source X-Ray 
Holographies: Bremsstrahlung, Fluorescence, and Multiple 

Energy X-Ray Holography 



Gerald A. Miller and Larry B. Sorensen 
Department of Physics, University of Washington, Seattle, Washington 98195 

Abstract 



Quantum electrodynamics (qed) is used to derive the differential cross 
sections measured in the three new experimental internal source ensemble x- 
ray holographies: bremsstrahlung (BXH), fluorescence (XFH), and multiple- 
energy (MEXH) x-ray holography. The polarization dependence of the BXH 
cross section is also obtained. For BXH, we study analytically and numeri- 
cally the possible effects of the virtual photons and electrons which enter qed 
calculations in summing over the intermediate states. For the low photon and 
electron energies used in the current experiments, we show that the virtual 
intermediate states produce only very small effects. This is because the un- 
certainty principle limits the distance that the virtual particles can propagate 
to be much shorter than the separation between the regions of high electron 
density in the adjacent atoms. We also find that using the asymptotic form 
of the scattering wave function causes about a 5-10% error for near forward 
scattering. 
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I. INTRODUCTION 



Fifty years ago, Gabor proposed electron holography as a method to improve the res- 
olution of electron microscopes so that atoms could be directly imaged ffl. Gabor 's idea 
was to focus the electron beam to a very small region of space just outside the sample to 
produce a nearly point source of radiation, and to record the interference pattern between 
the spherical reference wave from this point source and the spherical object waves produced 
when the electrons scattered from the atoms in the sample. This photographically recorded 
interference pattern would then be used as the diffraction grating in an optical reconstruc- 
tion system. Although Gabor's dream to directly image atoms using holographic electron 
microscopy has never been realized (because the quality of the best electron lenses is only 
about as good as that of a raindrop for visible light ||), Gabor's suggestion produced the 
optical holography revolution with the advent of lasers to provide the necessary coherent 
monochromatic external reference waves. 

Ten years ago, Szoke pointed out that the necessary coherent spherical reference wave 
could also be created by generating the electron reference wave inside the sample ||. In this 
case, the spatial coherence comes from the small spatial extent of the internal electron source. 
Szoke's internal source electron holography suggestion generated a flurry of activity 0], and, 
in the past five years, Gabor's dream of directly imaging atoms with electrons has been 
realized using photoelectrons || , Auger electrons |§ , diffusely scattered low energy electrons 
@, and diffusely scattered Kikuchi electrons §|. However, because electrons interact very 
strongly with atoms, the scattered object waves are not very good spherical waves (there is 
a strong angular variation of the magnitude and the phase of the electron-atom scattering 
amplitude), and multiple scattering produces "electron focussing" effects along the lines of 
atoms in the sample which are important. Consequently, these new electron holographies 
produce pictures of where the atoms are, but they do not accurately reconstruct the atomic 
positions. Szoke also proposed internal source x-ray holography using fluorescence x-rays. 
Because x-rays interact weakly with atoms, internal source x-ray holograms should produce 
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much more accurate atomic resolution images than the internal source electron holographies 
PJT0|j. Unfortunately, the price for this is that the modulation of the intensity in x-ray 



holograms (1—5x10 3 ) is about 100 times weaker than the modulation in electron holograms 
(1-5 x 10- 1 ). 

Earlier this year, the first atomic resolution x-ray holograms were produced using x-ray 
fluorescence holography (XFH) JTT|] and multiple energy x-ray holography (MEXH) ||12|| . 
About three years ago, stimulated by the advantages of multiple energy electron holography 
13| and the promise of XFH ||, we started developing a new kind of internal source x- 



ray holography which uses bremsstrahlung photons created inside the sample ||14|| . The 
primary motivation for this paper is to provide the theoretical foundation for experimental 
bremsstrahlung x-ray holography (BXH) starting from quantum electrodynamics. The bulk 
of this paper is devoted to BXH, but we also show how the same quantum electrodynamic 
foundation applies to XFH and MEXH. 

Bremsstrahlung x-ray holography is very attractive for three reasons: (1) Bremsstrahlung 
allows hard x-rays to be produced from low Z atoms. High quality holograms require the 
wavelength to be much smaller than the spacing between the atoms. The characteristic 
fluorescence energies of many interesting and important low Z elements are too low to 
provide good images using XFH. (2) Bremsstrahlung produces x-rays with a wide spread in 
their energy and allows multiple energy holograms to be recorded simultaneously by energy 
analyzing the bremsstrahlung photons. To accurately reconstruct a three-dimensional object 
in real space, we need information over a three-dimensional volume in reciprocal space. To 
overcome the problems in the internal source single-energy electron holographies, several 
multiple energy methods have been developed |13j. These multiple-energy electron methods 
eliminate the twin images, greatly reduce the effects due to the strong angular variation of 
the magnitude and phase of the electron-atom scattering amplitude, and reduce the noise 
in the reconstructions. The MEXH method was developed in analogy to the multiple energy 
electron methods to provide higher quality holograms than the single-energy XFH method. 
(3) The bremsstrahlung production cross section is extremely high. A conventional 400 watt 



x-ray source produces about 10 13 short wavelength bremsstrahlung photons per second into 
47r steradians ||14|| . If all of these photons could be collected and energy analyzed, a very 



high quality BXH could be generated with a tabletop apparatus in a few hours. 

The implementation of bremsstrahlung x-ray holography raises a number of interesting 
theoretical questions. We first recall that bremsstrahlung photons can have any energy from 
nearly zero to the energy of the incident electron, and the spectral intensity diverges at 
low energies: I{uo)duj ~ uj^ 1 duj. The bremsstrahlung photons, which scatter in the target 
crystal to produce the object waves, are intermediate or virtual particles. The momentum 
of each intermediate bremsstrahlung photon can take on any value; we must integrate over 
all virtual momenta in computing the scattering amplitude. Of course, this must occur 
for any intermediate particle that produces an object wave, but the possible problems are 
potentially more serious here for the virtual bremsstrahlung photons because of the broad 
nature of the bremsstrahlung spectrum. 

More generally we should ask: Are quantum mechanical effects ever important in internal 
source x-ray holography? Or does the simple wave picture always work? If a quantum me- 
chanical approach is needed, what is the correct quantum mechanical description of internal 
source x-ray holography? When can multi-path photon interference be treated by the scalar 
wave equation approximation to Maxwell's equations instead of the full theory of quantum 
electrodynamics? To answer these questions, we develop a quantum electrodynamic treat- 
ment of the three internal source x-ray holographies, BXH, XFH and MEXH, and compare 
it in detail with the simple wave picture. 

It is useful to provide this connection between the fundamental theory (QED) and these 
new holographies. Almost all of the work in this field has been based on the simple wave 
picture. For example, Barton's original holographic inversion procedure [|15j for electrons 
is based on the Helmholtz-Kirchoff inversion procedure for classical scalar waves. However, 
in the non-relativistic limit, Rous and Rubin [IBJ have recently shown how the Lippmann- 
Schwinger equation can be used to provide solutions to the Schrodinger equation which 
correctly describe the physics of the single-energy electron internal source holographies. 



For BXH in particular, classical electrodynamics will not produce the correct answer 
at high energies because the intermediate photons and electrons are virtual: the square of 
their four momenta may not be equal to the square of their rest masses, p 2 7^ m 2 . The 
physics can be divided into on-the-mass-shell amplitudes (called the "on-shell" or "real" 
processes) when p 2 = m 2 , and off-the-mass-shell amplitudes (called "off-shell" or "virtual" 
processes) when p 2 7^ m 2 . We show by explicit calculation that the virtual photons and 
virtual electrons do not propagate over the entire distance between the regions of high 
electron density in two adjacent atoms, and consequently classical electrodynamics predicts 
the right behavior. The reason for this comes from the uncertainty principle. The amount 
of off-shellness, or virtuality, of photons of energy k Q and momentum k is measured by the 
quantity k% — k 2 = Q 2 . If k > \k\, the virtual photon is not massless and its range is 
Q^ 1 which is much smaller than the interatomic spacing a. If k Q < \k\, Q 2 is negative and 
exp(i\Q\r) oscillates rapidly for r ~ a and any important contributions are cancelled. 

Thus our main result is that for real atoms in real solids, excited to emit bremsstrahlung 
or fluorescence x-ray radiation, classical electrodynamics works very well because there is 
no significant overlap between the regions of high electron density in the adjacent atoms. 
However, because the intermediate state photons and electrons in the internal source x- 
ray holographies are virtual, it is important to use quantum electrodynamics to derive the 
equations necessary to analyze the holograms. We provide that derivation for BXH, XFH 
and MEXH. 

The separated atom approximation, that we use to show that the full quantum electro- 
dynamic treatment reduces to the classical electrodynamic expressions for real atoms in real 
solids, is formally analogous to the separated scatterer approximation used in analyzing high 



energy hadron-nuclear scattering experiments p~7|JT8| JTl 

Almost all of our knowledge of the atomic scale structure of bulk condensed matter has 
been determined from measurements of the quantum mechanical interference patterns that 
arise from particle-crystal scattering. How do the new x-ray holographies compare with 
crystallography, and what are the other possibilities? There are four equivalence classes of 



quantum mechanical interference patterns that have been used to determine structure: (1) 
In crystallography, there is an external source of particles which are sent into the crystal in 
nearly plane wave states. In the usual kinematic scattering limit, these particles coherently 
single scatter from many atoms in the crystal. The interference between these many single- 
scattering events produces the Bragg peaks |2TJ. (2) In internal source holography, there is 
an internal source of particles which leave the crystal in nearly spherical wave states. These 
particles coherently single scatter from the object atoms in the crystal. The interference 
between each of these single-scattering events and the strong direct path reference beam 
produces the Gabor zone plates in the hologram. (3) In external source holography, there 
is a coherent external source of particles which is sent into the crystal in nearly plane wave 
states. These particles coherently single scatter from the object atoms in the crystal. The 
interference between each of these single-scattering events and the strong reference beam 
produces the hologram. Unfortunately, the necessary coherent hard x-ray sources are not 
yet available [^2| and when they become available they will probably destroy the sample in 
the process of making the hologram [53]. (4) In the Kichuchi and Kossel methods, there is 



an internal source of particles which leave the crystal in nearly spherical wave states. These 
particles coherently multiply scatter from many atoms in the crystal. The interference 
between these many multiple-scattering events produces the Kichuchi and Kossel patterns. 
These multiple-scattering patterns also contain useful holographic-like information p4j , but 
this information is different than the single-scattering holograms. 

The remainder of this paper is organized as follows. Sect. II outlines the standard 
classical scalar wave derivation of the intensity of the holographic interference pattern for 
the interference between a spherical reference wave and a spherical object wave. Sect. Ill 
is devoted to deriving an expression for the corresponding cross section for the intensity 
of the holographic interference pattern for bremsstrahlung holography. Since this paper is 
concerned with the possible effects of virtual photons, it is sufficient to consider the process 
as being bremsstrahlung production by the source atom followed by photon scattering by 
the object atom, located at a displacement r from the source atom. We find that we can 



simplify the expression for this cross section and apply it to holography if the atoms can be 
treated as well separated so that only real photons propagate from the source atom to the 
object atom. We then show for real atoms in real crystals that the regions of high electron 
density are sufficiently well separated. Our separated atom approximation is presented in 
Sect. IV. The bremsstrahlung energies for the experiments we are considering are 40-60 
keV and at these energies, the x-ray-atom scattering cross section is dominated by the 
Thompson process. So we study the photon virtuality effects for bremstrahlung production 
followed by Thomson scattering first. For this case, the corrections to our separated atom 
approximation are defined and shown to be entirely negligible in Sect. V. Near resonance, 
the anomalous scattering amplitude can become comparable to the Thompson scattering 
amplitude, and these two amplitudes interfere. We consider this case in Sect. VI, where 
we use the numerical results of the previous sections to justify the immediate use of the 
separated atom approximation. In XFH, the photon is produced by fluorescence radiation, 
where the excited atomic state is produced by electron or photon impact, and the atom 
decays via photon emission. The emitted photon can be scattered by another atom to 
produce an object amplitude which will interfere with the direct reference amplitude. This 
is discussed in Sect. VII, where the necessary amplitudes for this process are presented. In 
MEXH, a real photon is sent into the sample from outside. This photon has a direct path 
to the detector atom and a collection of single scattering paths to the detector atom via the 
object atoms which will interfere with the direct path. This is discussed in Sect. VIII. The 
final Section is devoted to a brief summary and discussion. 



II. CLASSICAL INTERNAL SOURCE ENSEMBLE X-RAY HOLOGRAPHY 

When an x-ray photon is created inside a solid and is detected outside, the quantum 
mechanical interference between the different paths that the photon takes as it leaves the 
solid will produce a holographic image of the atoms around the position where the photon is 
created. As we show below, the probability distribution for the photon intensity is a Gabor 
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hologram. In contrast to the usual external source x-ray holography where the reference wave 
comes from outside the sample, in internal source x-ray holography, the wave corresponding 
to the direct amplitude (i.e., the amplitude for the photon to leave the solid without any 
interactions) serves as the reference wave, and the wave corresponding to the amplitude 
produced by single photon-atom scattering plays the role of the object wave. 

Because the amplitude for photon-atom scattering is weak for hard x-rays, the reference 
wave will be much stronger than the object waves, and this strong reference wave limit is 
the ideal holographic situation because the hologram is then dominated by the interference 
between the reference wave and the singly scattered object waves. In this limit, the inter- 
ference between one object wave and another object wave, and the interference between the 
reference wave and the low order (double, triple, ...) multiple scattering object waves is 
much weaker than the interference between the reference wave and the single-scattering ob- 
ject waves. In nearly perfect crystals, the interference between the high order object waves 
can become comparable to the reference wave and is responsible for the Kossel (x-ray) and 
Kikuchi (electron) patterns. However, these multiple scattering features are sharp in angle 
and therefore can be easily removed from the hologram. 

To develop the simple classical wave picture for internal source ensemble x-ray holog- 
raphy, consider first just the two atoms shown in Fig. 1. The full internal source x-ray 
hologram can be obtained from this two atom case by summing over all object atoms for 
each source atom, and by summing over all source atoms. Suppose the source atom at the 
origin emits radiation which is detected in the far field, and suppose also that prior to de- 
tection the radiation is scattered by a second object atom located at position a. The direct 
and single-scattering paths produce an interference pattern. When the polarization is not 
important, this problem can be treated as due to the interference between two scalar wave 
fields, with the scalar field representing a component of E or B. 

In this approximation, the first atom emits a scalar spherical reference wave R of the 
form 
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and the second atom emits a scalar spherical object wave O of the form 

pika pik\r— a\ 

o = — f a (e) e —— , (2) 

a \r — a\ 

where again a is the position of the second atom. In the far field, r ^> a, the composite 
amplitude M = R + O takes the form 

M = — (1 + e ika e~ &s ) (3) 

r a 

where f a (Q) is the standard atomic scattering amplitude for an incident plane wave, and 
we have followed the tradition in this field of neglecting the higher-order spherical wave 
corrections. 

The far field intensity is given by the square of the composite amplitude M 

I = M*M = R*R + R*0 + RO* + 0*0. (4) 

Holography records much more phase information than crystallography, but there is still a 
"holographic phase problem" present in the holographies that can only measure the inten- 
sity (e.g., laser, electron and x-ray holography) and not the amplitude (e.g., acoustic and 
microwave holography). The intensity holographies record both the information we want 
about the object in the R*0 term, and a copy of the complex conjugate of this information in 
the RO* term, which produces a non-existent twin to the object during the reconstruction. 

The far field object plus twin holographic interference cross term R*0 + RO* is propor- 
tional to 

Re [f a (Q)] cos(ka — k ■ a) — Im [f a (Q)] sin(ka — k ■ a). (5) 

In this paper, we treat the case in which the source atom emits x-rays with wavelength A. 
Note that to obtain significant holographic oscillations, it is important to make A < a, and 
that the best holograms will be produced when A C a. 
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Although all of the existing and proposed internal source holographic techniques actually 
depend on the multi-path quantum mechanical interference of the particle emitted by the 
sample, the essential features of the holograms can be (and have been before this paper) 
obtained from the simple wave picture of the process outlined above. 

There are two essential ingredients of atomic resolution internal source holography: (1) 
There is a localized source inside the sample. This localization provides the necessary spatial 
coherence of the source. The particle can be localized by being created inside the sample— 
this is the case for the bremsstrahlung and fluorescence x-ray holographies described in 
this paper. The particle can also be localized by being ejected from a specific quantum 
state in the sample — this is the case for photoelectron and Auger electron holography. In 
addition, the particle can be localized by an incoherent inelastic scattering event — this is the 
case for diffuse LEED holography and diffuse Kikuchi electron holography. The analogous 
incoherent scattering localization is possible theoretically for photons via thermal diffuse 
x-ray scattering and via Compton scattering. (2) There is interference between the direct 
reference wave and the singly scattered object waves. This requires a coherent scattering 
event at the object atoms. If the object atoms scatter incoherently, the interference in the 
final state will not occur. 

We shall show how these features arise from quantum electrodynamics. 

III. BREMSSTRAHLUNG X-RAY HOLOGRAPHY (BXH) 

In the bremsstrahlung process, an electron incident on a solid radiates a photon: e(pi) + 
solid — > e'ipj) + j(k) + solid', excited. The effects of thermal motion are ignored. In this 
paper, we consider the case where the initial and final electronic states of the atoms are the 
same, and where the solid is a collection of fixed atoms, i.e., we do not consider the effects of 
thermal motion. The Feynman diagrams for bremsstrahlung holography are shown in Figs. 2 
and 3. The complete holographic amplitude M. is the sum of the reference wave amplitude 
1Z, given by the crossed and uncrossed Born bremsstrahlung terms shown in Fig. 2, and the 
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object wave amplitude O, given by the crossed and uncrossed Compton scattering terms 
shown in Fig. 3: 

M=TZ + 0. (6) 

Here the QED reference TZ, object O, and hologram Ai amplitudes are analogous to the 
corresponding classical R, O, and M terms in Eqs. (P through ([5]). 

Before starting, it is useful to sketch the notational conventions used in this paper. The 
various four-momentum vectors are represented with Roman typeface, e.g., k, p«, p/, and 
the three-vector spatial components are indicated with italic typeface with explicit arrows, 
e.g., k, pi, and pf. The magnitudes of three-vectors are written explicitly, e.g., \k\. Thus 
for the bremsstrahlung reference amplitude shown in Fig. 2, the conservation of energy and 
momentum is written as p^ + q = p/ + k, where u = k°, p° = Ei, and p° = Ef. Here 
q is the four-momentum supplied by the target, q = (0,q), where the arises from our 
condition that no atoms (or nuclei) be excited. Three-momentum conservation is expressed 
as q + Pi — Pf + k. The notation and conventions of Bjorken and Drell are used 
throughout the present work, and our units are such that both h and c are unity. 

The cross section for the holographic interference pattern is related to the square of the 
holographic amplitude M. by: 



— = m 2 £1 2tt — — 6(Ei — m — u)) (7) 

dn dQ f d\k\ Pi 2 (2tt) 6 V ; V ; 

where m is the mass of the electron, Q represents the outgoing angles of the photon and 

Qf those of the electron (pf). The quantity \Ai 2 \ is obtained from \Ai\ 2 by squaring the 

magnitude of 1Z + O, summing over the spins of the final electron, and averaging over the 

spins of the initial electron. 

The reference amplitude R shown diagrammatically in Fig. 2 is evaluated as 



TZ{k,q) 



Ze p e 2 " < f ■ •■ ^ 



12 



mm w 



2pf ■ k 2pi ■ k J 2pf ■ k 2pi ■ k 
where e 2 /47r = a, the proton charge is the negative of the electron charge e p = — e, and 
F(\q |) is the electronic contribution to the atomic form factor, normalized so that F(0) = 1. 
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This matrix element is proportional to e M = (0, e). The nuclear form factor is essentially 
unity for the kinematic range of the current x-ray holography experiments. The deviation 
of the term (1 — F(\q |)) from unity represents the screening effect of the atomic electrons. 
If q = 0, the atom acts as a neutral object and there is no bremsstrahlung. 

The only significant approximation made in obtaining Eq. (Q) is that the initial and final 
state electron-nuclear Coulomb interactions have been neglected. The influence of these 
interactions, which can increase the value of the computed cross sections significantly, can 
be reasonably well approximated by multiplying the above amplitude by the product of 
the continuum electronic wave functions evaluated at the nuclear center — this is the Elw- 



ert approximation of Ref. [26[. This is a well motivated approximation classically because 
the acceleration that leads to the bremsstrahlung takes place in the vicinity of the nucleus. 
Detailed numerical studies |27| have confirmed the qualitative accuracy of the Elwart ap- 



proximation. Multiplying our amplitude by this factor does not influence the propagation 
of the virtual photon between atoms, which is our principle concern. Thus we shall ignore 
the initial and final state interactions here in our study of the potential off-shell effects. 

It is convenient to define the expressions in the bracket as e • B(k), so that the reference 
amplitude can be rewritten as 

Ze e 2 

1l(k,q) = -^(l-F(\q\))B»(ky. (9) 

As shown in Fig. 3, the virtual photon (k') is produced by the source atom, propagates to 
the object atom, located at a separation f from the source, which scatters the virtual photon 
k' so that the final photon k is produced. The object atom scattering is dominated by the 
Compton scattering of the photon by the atomic electrons. This is because the photon-atom 
scattering is larger than the photon-nuclear scattering, by the ratio of the proton mass to 
the electron mass for the Thompson term, or by the ratio of the squares of the atomic and 
nuclear radii for the dipole terms. The virtual bremsstrahlung matrix element is denoted as 
B^k') and the Compton rescattering transition matrix as C M (e, k, k'). The evaluation of the 
Feynman graphs shown in Fig. 3 uses standard techniques |25| . Here we also carry out the 
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integration over the time component of k', which gives us a delta function setting k f0 = u. 
Then we arrive at the expression for the object amplitude: 

rM. f , cvu v y\ » d3k ' ^ - F ^ 



O = -Z 2 e\ / < C^e, k, k') > ^ e^ k ~ k ^ /0 v ^ 2 177 (10) 

where q' = pf — pi + k'. Note that fc /2 = uo 2 — k ' 2 ^ 0. 
The bremsstrahlung matrix element is given by 

where standard spinor notation is used. Note that the Compton term < C M (e, k, k') > is the 
atomic expectation value of the virtual-to-real Compton transition matrix which converts 
the virtual photon k ' to a real one with four momentum and polarization (k, e) . Thus 

C{e, k, k') = ^S F (Pi -k)i+i S F (P f + k)Y, (12) 

where S^(P) is the relevant propagator for the bound electrons. For example, when the 
electrons are treated as free, the Compton transition matrix element is given by 

*'> s °> { Y r^j-irJ + ' u ' (13) 

where the upper case spinors Ui, Uf represent the initial and final states of the free electron, 
with Pi + k' = P f + k. 

But the atomic bound states are more interesting. We may better understand the oper- 
ator of Eq. ( |12"D by noting that in the relativistic theory, the origin of the Thompson term 
comes from the terms involving the creation of virtual electron-anti-electron pairs |28] . The 
resulting two-electron plus anti-electron states live only for a very short time, so that we 
may ignore interactions with the other particles of the atom. The remaining terms can be 
seen in the non-relativistic limit up to 0(p 2 /m 2 ) as arising from the two interactions of the 



dipole operator ||28|| . Then, we may write the Compton transition matrix element as a sum 
of terms so that 

< tf"(e, k, k') >=< T»(e, k, £;')> + < R"(e, k, k') > (14) 
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where the Thompson scattering is denoted as T M (e, k, k') and the resonant scattering and 
other contributions are denoted as < -R M (e, k, k') >. 

When the long wavelength approximation is valid, we obtain 

t^u, , 9 9 < *|e • D\n >< n\Dr\i > . 

< W{e, k, k) >= uj 2 e 2 j: ' ' p ,\" <W + 

n £ji + U) — hi n + n 

2 2 < >< n\e ■ D\i > . . 

we 2- p — - — p ^' ( 15 ) 

n Hii — LU — Hi n 

where the dipole operator D is given by 

z 



D = $>1 (16) 



i=i 



and s*j is the displacement of the i-th electron from the atomic center. The vector D is 
simply the sum of electronic dipole operators. The quantities E n and T n are the energy 
and the width of the excited state n. The result shown in Eq. ( |T3j) indicates that only the 
three-vector R part of R 11 enters into the expression for the amplitude. This is because 
each atomic photon emission/absorption is controlled by an e • D operator. For unpolarized 
atoms, R must be proportional to e, so that it is convenient to define a strength function 
S(u) such that: 

R = iS(u). (17) 



We use this definition along with Eq. (|i~5|) to obtain 
S(u) = u 2 e 2 \ < i\e ■ D\n > '° ' 1 



+ 



(18) 



Ei + u — E n + \iT n Ei — u — E n _ 

As noted in the introduction, we shall proceed by first studying the effects of Thompson 
scattering by the object atoms. An explicit evaluation yields 

<r>=-^re ! F(|l-fc'|), (19) 

where for simplicity we take the scattering object atom to be of the same type as the source 
atom which produced the virtual photon. In this and the following two sections we shall 
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consider photons for which the Thompson term is dominant. We shall return to the resonant 
corrections to the Thompson term in Sect. V. 



We compute \M. 2 \ by squaring IZ+O, keeping only the Thompson scattering contribution 
in O, by summing over the final electron spin, and by averaging over the initial electron 
spin. The result is 



1 



l*»l = 5 E 

SfSi 



fZe 3 . \ 2 Z 3 e 8 



d*k> (1-F(\q'\) , ( ^ iyP F(\k-k>\) 



x 2 Rel B*(k) ■ e / -— - ^ — ^ u e 



(2vr) 3 |gf CJ 2_ A ./ 2+ie 
x < T M (e, fc, fc') > 5 M (/c')|J . (20) 

Here the term of order Z 4 e 10 , which is much smaller than the lower order terms, has been 
neglected. 

The result given by Eq. (^), as specified by the matrix elements given by Eqs. (|TT|) and 
(|19|), is our main result. It gives the bremsstrahlung holography cross section when the ob- 
ject scattering is dominated by the Thompson term. We shall make a further simplification 
to facilitate a first evaluation: we will keep only the numerically most significant terms of 
B and B, i.e., the ones proportional to e • p) and e • p"j. This leads to the standard classi- 



cal expression for the bremsstrahlung cross section ||29|| . We have explicitly evaluated the 
neglected terms numerically and found that their neglect produces an approximately con- 
stant 10% reduction in the cross section over the kinematic region where the bremsstrahlung 
holography experiments will operate. Our numerical results can be understood by noting 
that $7° ~ io ■ e x k. This spin-dependent interaction is a magnetic effect proportional to 
V x A, which therefore does not interfere with the terms we keep. Furthermore, the two 
spin-dependent terms of Eq. (§) partially cancel and their sum is smaller than the leading 
term by about k/m. 

We carry out the average over electron initial spin, and the sum over final electron spin. 
The result for the bremsstrahlung holography cross section is: 

d3(X .2 b/l 1 w _ ATT7T9 1 roPP i_ _ , i„2i 



m -yirr - 9(Ei -m-iu)\M\ 2 — [8EiE f - A Pi ■ p f + Am 2 ] (21) 



dttdtt f d\k\ \Pi\ 2 (2tt) 5 
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m 2 



where 



\M\ 2 ={^(l-F(\p f -p i + k\))) [e ■ V(k) e ■ V(k)}^ 2(1 F{Wl)) Re I(k, r) . (22) 
\ \Q I / 19 1 m 

Here the quantity I(k,r) is given by 



, r, m .-nf-ty? < V - n\f'\)) F(\k-kf 



m r> m ,. m j^i. vw ^F7t} • (23) 

where q' = Pf—Pi + k' and q = pf — pi + k. The vectors V(A;) and V^A;') are given by 

^(fc) = JL-JL (24) 

P/ ■ K Pi • K 

= 2p f ■ k' + h* ~ 2 Pl ■ k>- Jfc« ' (25) 



The cross section for the intensity of a bremsstrahlung hologram given by Eq. fl2~T|) 
together with the definitions given by Eqs. (]22"D through (|25| ) is the complete solution to the 
bremsstrahlung holography problem. What remains to be done is to carefully analyze these 
equations to see how the classical holography equations emerge in the classical limit, and to 
see how large the quantum effects are, and when they are important. That is the content of 
the next three sections. 



IV. SEPARATED ATOM APPROXIMATION 

The goal of the bremsstrahlung holography experiments is to determine precise infor- 
mation about the location of the object atoms, which is represented in Eqs. ( |2T| ) through 
by the vector r. The standard holography expressions involve an interference term of 
the general form j3 e tujr e - %k - r j r where (3 is a known function. A quick look at Eqs. ( |2"2"| ) and 
(p3|) could lead one to dismay. How could that integral ever have the simple form required 
for holographic investigations? We indicate a solution by considering the situation when 
the two atoms are very far apart, i.e., in the limit where r approaches infinity. Here our 
intuition provides a guide: the process must proceed by bremsstrahlung from the source 
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atom followed by photon propagation along the direction of r and Thompson scattering by 
the object atom. Thus the bremsstrahlung makes a real photon with momentum k = tor 
and energy u, and the Thompson scattering changes the direction of k to k. 

In this case, the photon has four momentum k = (u, k) and k ■ k — 0. The propagating 
photon is on-shell. This situation is simple, but the full integral of Eq. (^) is not. However, 
we will evaluate this integral by developing expansions in which the leading term is correct 
in the limit that r is very large. We shall keep the leading term and the most important 
corrections. We call this approach the separated atom approximation. In practice, this 
amounts to replacing k ' by n in certain terms in the integrand. For example, we shall show 
below that replacing V(k ') by V(k) and q' by pf — p^ + k are excellent approximations. 

We may use the uncertainty principle to better understand why the propagating photon 
must be real for infinite values of r. Our virtual photons have energy u, but the magnitude of 
the three momentum varies from to infinity in the integration. Let us define Q 2 = uo 2 — \ k '\ 2 
to provide a measure of the violation of energy conservation required to make the virtual 
photon. This is simply the square of the energy-momentum four vector, which vanishes for 
real photons. If Q 2 < the wave is a decaying exponential of the form e~' a< ^ r '/r, which has 
a small value. The interpretation of this Yukawa form is that the photon lives for a time 
h/Q, so that its maximum range is hc/Q. If Q 2 > 0, the wave is of the form e in ® r /r. The 
effect of this term is very small because of the oscillations of the integrand in the integral 
over d 3 k '. Thus, the net result is that only the real photons with Q = reach the object 
atom. 

Are the atoms in a real solid sufficiently separated so that all the virtual photon effects 
are gone before the bremsstrahlung photon reaches the nearest atoms? We argue that the 
answer is yes, at least for most solids at typical experimental bremsstrahlung holography 
energies, by considering the specific example of crystalline copper. In crystalline copper the 
atoms are separated by a distance Rn ~ 2.5561 A w 4.83ao where ao is the Bohr radius 
~ 0.529 A. Where are the electrons in each atom? The electron density for isolated copper 



atoms and for crystalline copper calculated using the FEFF computer code |30| is shown 
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in Fig. 4. Note that this density is sharply peaked at small values of s since most of the 
electrons are within 1 A of the nuclear center of the atom, and that the electron densities 
for isolated atoms and for atoms embedded in the solid are very similar. In particular, the 
root mean square radius of the displayed density is 1.08 do- Thus the closest separation r 
between the copper atoms is about 5 times the typical value of the distance s between an 
electron and the nucleus. At the very least, it is reasonable to expect that the separated 
atom approximation is a good starting point. 

Our procedure is to examine a set of approximations to the full results for I(k, r) given 
by Eq. (0). We will define the on-shell, separated atom approximation as the result of 
setting k' = k in Vi(k') and in q'. Then the on-shell approximation I on (k, r) to the full on- 
and off-shell integral I(k, r) is given by 

, = (l-F(P /: a + g|))(l-F(|f|)) ■ . y (k ,.y (K) £ ^ (26) 

\Pf ~Pi + K\ 2 

where 

J {^) u 2 -k' 2 + ie 

Again, the subscript "on" is to remind us that setting k' equal to k causes the propagating 
photon to be on-shell. Its virtuality has decayed by the time it reaches the near-neighbor 
atoms, and the square of its four momentum has vanished. The next section (Sect. V) is 
devoted to the demonstration that the on-shell I on (k,r) given by Eq. (|26|) is an excellent 
approximation to the full I(k,f) given by Eq. (|23|) . 

The first step is to understand the integral J on . We can gain some insight by converting 
the integral over the momentum into one involving positions. We use 

F{\k-k'\) = J d 3 s p(s)e-^-^- y , (28) 

so that the on-shell J on integral given in Eq. (pTD simplifies to 

J m {k, r) = -— / d 3 s p(s) ——e-^ s+ ^. (29) 
47T J \s + r\ 



If r ^> s for the important regions of p(s), we may replace e l0J ' s+ ^/|s + r\ by e luir e luir ' s /r. 
This gives 

lim J on (A;,f) = F(\k - k\), (30) 

) — i-oo 47T r 

which has the usual classical holographic form. The spherical Green's function e tujr /r cor- 
responds to the form of the wave at infinity. This form arises from the pole in the integral 
for J on at \k'\ — u). 

To check the asymptotic approximation given in Eq. (^0|) we numerically compared the 
full J on given by Eq. ( p9[) with its approximation given by Eq. (|30|). The results are shown 



in Fig. 5. Note that the approximation is excellent except when k \\ f. Even then, the 
full theory and the approximation produce very similar holographic interference patterns; 
the asymptotic approximation produces a pattern about 10 % smaller than the full theory 
when k \\ r. This agreement between the full quantum electrodynamic calculation and the 
simple classical holography equations shows that bremsstrahlung holography is possible: 
the reduction for k || r does not significantly change the oscillatory form with its strong 
dependence on k • f '. 

There is no need to use the asymptotic approximation in numerical work. We may use 
the correct value of J on and maintain the explicit holographic form. This involves expanding 
the form factor in terms of Legendre polynomials Pl(/c • k'): 

F(\k-k'\) = Y,F L {u J r)P L {k ■ k'), (31) 

L 

where uj = \k\. Combining the partial wave expansion for the atomic form factor given by 
Eq. (|3l|) with the full expression for J on given by Eq. (|27D yields the partial wave expression 
for J n 



' on 



J on (k,r) = iuJ2i L FL(u;r)h { l\L;r)P L Ck ■ r), (32) 



L 

,(1) 



where h L (ur) are the outgoing spherical Bessel functions. All that is required for this to 
hold is that r be bigger than the maximum value of s (i.e., circa 2.56 A for copper) occuring 
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in the integral (|29|) . Since r = 2.56 A and Fig. 4 shows that p(s) is less than 10~ 3 of its 
maximum value for s > 0.5 A, this condition is met. We may also understand the relation 
between this expression and its limiting form shown in Eq. (|30|) . The use of the asymptotic 
form of the outgoing spherical Bessel functions: 



Jim4 1) (x) = H) i+1 ^, (33) 



leads immediately to the result shown in Eq. (|30|). The corrections to this asymptotic form 
are thus of order 1/x times the original result. Thus we see that the expected first correction 
to Eq. (pO) is of the order of l/wr « 1/25 for uo = 20 keV. 



Equation ( |32| ) allows us to understand why the difference between the asymptotic 
approximation given by Eq. (|30f) and the exact result given by Eq. (|3~2"D is largest for 
k ■ f = 1. The terms FL{uor)hL{uor) monotonically approach zero as L increases. The 
function Pi{k • r — 1) = 1 for all L, so that the terms with large L (these are the terms for 
which the approximation generated by Eq. (|33|) is less accurate) add constructively. This 
can also be seen (without using the partial wave expansion) by examining the integrand of 
Eq. fl29|). If k ■ r = 1 {i.e., when k || f), the term uo\s + r\ — k ■ s is greatly reduced so that 
the contributions of the larger values of s are less inhibited by the oscillating exponential 
than for other values of k ■ r '. 

The partial wave expression (p2|) systematically gives all of the corrections to the classical 
holographic form given by Eq. ((3^). However, it is useful to provide another approximation 
which shows us why the relevant integrals are dominated by terms in which k = k. The 
idea is to approximate F(\k — k'\) using 

F{\k-k'\) fa F( qi ) + {%' - k) ■ V Ql F( qi ) (34) 



where q\ = \k — k\. We use (|34]) in the integral (|27|) and note that the k! appearing in the 
numerator of the integral (^) can be replaced by a gradient on r. Thus we find 

^ i p iuir 

Jon(k, r) « e~ ik * F( gi ) + 8J m , (35) 

47T r 
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with 

6J m (k, r) = -e- ltr V qi F{ qi ) ■ V on (k, f), (36) 

where 

™*'$-*)/&*'{^}- 

One may use the Legendre partial wave expansion of Eq. (|32"D above to obtain a more 
detailed expression for V on (k,r). But the main point is that the long distance behavior of 
the integral is that of e tujr . Since (— zV r — R)e luir = 0, the correction to the separated atom 
approximation must have an extra factor of order 1/ur ~ 1/25. We shall use expansions 
similar to that of Eq. (|34l) to systematically understand the short distance terms. The 
integral V on (k, r) will appear again. Furthermore, we shall often employ the technique of 
writing a complete expression as its on-shell approximation plus a term which is proportional 
to (—iV r — and vanishes in the asymptotic limit given by Eq. (|30|). 

V. SHORT RANGE TERMS 

— * 

In the previous section we showed how keeping the effects of the pole at \k'\ = uo led to 
the term with the long distance propagation. Here we show that this pole dominates the 
complete expression given by Eq. (|23|) . The vector k' appears in three places in this integral, 
in (1 — F(\q'\)), in l/\q'\ 2 , and in Vi(k'). We will denote these terms as the screening 
correction, the Coulomb photon propagation, and the electron propagation. We shall study 
each one separately. 



A. Screening correction 



Keeping the k! in the screening term leads to the integral, I s : 



J -(M= / 77T~T3 e | (1 ~ F(\q'\)) " 1 . (38) 
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Recall that q' = pf—pi + k'. Pole dominance of the integral would allow us to replace the 
k' appearing in q ' by k. Thus the integral may be approximated by using 

F(\q'\) = F(\p f - pi + k + {%' - k)\) « F(\p f -£. + £|) + (£'-£) . V g2 F(q 2 ), (39) 

where q 2 = \pf— Pi + k\. Using this in Eq. ( |3~8D leads to the appearance of k! in the numerator 
of the integral, which can again be replaced by a gradient on r. The result is 

I s (k, f) « (1 - F(q 2 )) J on (k, f) + 5I s (k, r) (40) 

with 

5I s (k, r) = -e~^V q2 F(q 2 ) ■ V on (k, r). (41) 

The leading long distance behavior of the integral of Eq. (|37|) required to evaluate V on 
is that of e lujr . But (— iV r — t£)e wr = 0. Furthermore, V q2 F(q 2 ) is of order (n/q 2 )F(q 2 ), 
with n ~ 4. Thus, the 5I s (k, r) screening correction term of Eq. (fiO| ) provides a correction 
which has an extra factor of F(q 2 )/q 2 r compared to the leading term. However, F(g 2 ) is very 
small, 1% at most, for the kinematics of this experiment. For typical kinematics q 2 ~ 100 
keV so (n/q 2 r)F{q 2 ) (4/125)(l/100) ~ 3 x 10~ 4 . The correction to the separated atom 
approximation due to the screening term is completely negligible here. This is shown in 
Fig. 6. 



B. Coulomb photon propagation 

If we keep the k' in the Coulomb photon propagator 1 / 1 q ' \ 2 , we need to evaluate the 
Coulomb integral I cou i(k,f): 

«*. - I ^ **** {^-^ ■ (42) 



(2tt) 3 I \p f + k'\ 2 u 2 - k' 2 + le • 

In this case, there are two sets of poles. One is the usual one at \k'\ — u, but there is 
also a set of poles off the real axis (in the complex \k'\ plane) corresponding to the zeros of 
\f>f — Pi + k'\ 2 . It is desirable to handle these pole terms separately, so we use the identity 
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1 1 

A ' B 



1 1 

A + B 



A + B 



(43) 



in the form 



\Pf -Pi + k 



'I 2 u 2 -k' 2 + ie 



+ 



k' 2 + ie \pf — f>i + k 



'12 



(44) 



u 2 + (pi - Pf) 2 + 2k' ■ (p f - p^ 
The first term has the pole at \k'\ = uo which is responsible for the long distance photon 
propagation. The second term has the above mentioned poles of \k'\ off the real axis. The 
vanishing of the denominator u 2 + (pi ~Pf) 2 + 2k' ■ (pf —Pi) occurs only when the two terms 
in the bracket cancel and causes no mathematical difficulty. The first term of Eq. fl4"3]) can 
be approximated by using k' = k, + [k' — k) and expanding so that 

1 1 11 



u 2 - k' 2 + ie u 2 + {pi - p f ) 2 + 2k' ■ (p f - p^ ~ u 2 - k' 2 + ie (Pf ~Pi + 

2(k'-K)-(p f -p l ) 
{Pf-Pi + k) 2 



(45) 



The first part of Eq. (45) corresponds to the separated atom approximation. The next 
term involves (Ar — k,) which yields the integral V on {k } f) of Eq. d37|) . The specific correction 
to J on is denoted as SJ^° ul which is obtained by keeping the second term of Eq. (55) in the 
I c integral given by Eq. ( ftS|) so that 



SJZ u \k,r) 



-2) 



d 3 k' c _itf_p w ( (k'-K)- (p f - p{) 

(2vr) 3 I ^ - fc/2 + ie 



(46) 



The k ' — k term of Eq. (H3) again leads to an extra factor of l/q^f, as compared to the 
leading term. The expected small size of this correction is confirmed by numerical evaluation. 
Indeed the second term of Eq. ( fl~5]) is negligible except for the angles for which the leading 
term vanishes. See Fig. 7, which shows the relative sizes of J on {k, r) / (pf — Pi + k) 2 and the 
correction to it due to 5J^ ul (k, f)/{pf — pi + R) 2 . 

What about the second term of Eq. (^4[) ? This is given by 



1 



1 



\pf -Pi + k'\ 2 \ u 2 + (pi - p f ) 2 + 2k' ■ (p f - p^ 
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It is necessary to treat the k! ■ k! terms correctly, but the term k ' ■ (pf —pi) may be evaluated 
by using k' — k + (k — k) and treating the difference term as a perturbation in k' — k. Thus 

1 ^ ^ 1 ^ 2(fc / -K)-(p / -pQ ^ 

|p) - Pi + fc'| 2 (pf - pi) 2 + 2k • (p f - p^ + k' ■ k' ((p f - p^ 2 + 2k ■ (p f - p^ + k' ■ k') 2 

If we treat this as a function of \k'\, the poles in the exact expression and in its approximation 
given by the first term, appear at the same positions. The second term can be thought of as 
correcting the value of the residue at the pole. Furthermore, it vanishes for well-separated 
atoms. Thus neglecting the second term is a good approximation. We use similar logic to 
write 

(48) 



(47) 



uj 2 + (pi - p f ) 2 + 2k' ■ (p f - ^) uj 2 + (pi - p f ) 2 + 2k ■ (p f - ^) 
We use the first terms of Eqs. ( |7]) and ( [48]) to estimate the second term of Eq. (S3) . We 
immediately expect that this second term is completely negligible because it has the form 
of the Fourier transform of 



k'-k' + P 2 

where P 2 = (pf — Pi) + 2k ■ (pf — pi). This Fourier transform falls off very rapidly with 
r, i.e., as exp(— Pr)/r. For typical values of P of about 100 keV ~ 50A , we will have 
Pr « 25, where a is the separation (circa 2.56 A). The exponential damping factor destroys 
this second term. This intuitive conclusion is also confirmed by numerical evaluation, but 
the strong exponential damping inherent in this term caused the effects of this correction to 
be too small to be plotted. 



C. Electron propagation 

The full expression for the Feynman graphs in Fig. 2 allows a new type of term, one in 
which the electron propagates over the distance r. The mathematical origins of this effect 
are in the poles of the electron propagator shown in Eq. (^) which arise via the appearance 
of the four-momentum of the virtual k ' 2 ^ in those denominators. There are two terms in 
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that equation, one arising from the uncrossed graph (Fig. 2a) and the other from the crossed 
graph (Fig. 2b). We shall study these terms in sequence, using our standard technique of 
writing k! = k+ (k' — K) and treating the second term as an expansion parameter — whenever 
possible without destroying the correct analytic structure. 



1. Uncrossed term 



Suppose we keep the full uncrossed term. This means that we must evaluate the integral 
Ji(k,r): 

d 3 k f e -i(k-k')-r 1 F(\k-k'\) ^ 



Jl(A; ' r ) ' (2tt) 3 2p f ■ k' + k' 2 + it k' 2 + it - k' 2 + it 
The product of propagators can be written 

1 1 



2p f ■ k' + k' 2 k' 2 + 



it 



2p f ■ k' 



(50) 



_k' 2 + it 2pf ■ k' + k' 2 _ 

The l/(k' 2 + it) is the photon propagator and we denote its contribution as the photon 
propagation term; similarly, the second part is the electron propagation term. The last 
factor l/2pf ■ k! vanishes only when the term in the bracket vanishes and so causes no 
mathematical difficulty We may then study two separate integrals 



J 1 (k,r) = K 1 (k,f)+K 2 (k,r), 



(51) 



where 



d 3 k' e -<(*-*V 1 Ffljfe-jfe'l) 



(2tt) 3 2p r k> k' 2 + tt t0 2-k>2 + l 



(52) 



it 



and 



K, 



d 3 k' 



-i{k—k')-r 



1 F(\k — k'\ 



(2tt) 3 2p f ■ k' + k' 2 + it 2p f ■ k' uj 2 - k' 2 + it 
We work first with K\. The manipulations are simplified by using 



(53) 
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1 1 , A 1 

+ ;2j9/ • (k - fe ) 



(54) 



2pj • fe' 2pj • k ' 2pf • k r 1 v 2pj ■ k' 

where the four vector k = (uj,k). Using this relation in Eq. (|52|) enables us to derive a 
differential equation for K\\ 



ATi(fc,f 



-2p f ■ 



k — — — k 



#i(Jfe,r 



(55) 



2pf ■ K 2pf ■ K 

We see that the first term is the separated atom approximation for this particular term. 
The quantity K, + iV r — k vanishes when acting on e r e ia ""/r, so that the second term 
is a correction. The effect of this term can be estimated by replacing K\ on the right hand 
side by J on /{2pj ■ k). Thus 



(2 Pf ■ k) K x (k,r) « Jon{k,f) + 6JZ cr (k,r), 



(56) 



with 



-2p f 



k + iVr — k 



Jon(k, r ) . 



(57) 



A brief calculation shows that once again the correction is proportional to the vector 
integral V on (k, r) of Eq. (j37D , and is down by about 1/ur compared to the first term. Explicit 
numerical evaluation confirms this estimate, the correction term J^ cr is indeed negligible, 
as shown in Fig. 8. For general purposes, it is useful to note that Eq. fl55|) has the formal 
solution 



ATi(fc,f 



1 



(58) 



2pj ■ k + 2pf ■ (k — k + i¥ r 
This formal solution gives us a controlled way to study some of the corrections to the 
separated atom approximation. 

The K2 term given by Eq. ( |53| ) represents new physics occuring in this two-atom process. 
To see this, recall that 2p f ■ k' + k' 2 + ie = (p f + k') 2 - M 2 + ie = (E f + uj) 2 - m 2 - 
(pf + k') ■ (pf + k') + ie. We write this in terms of a four vector W=(Ef + u,pf + k') as 
2pfk' + k' 2 + ie = W 2 — m 2 +ie. The zero in this term is a pole in the integrand representing 
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the long distance propagation of the electron of four- momentum Won = (Ef + uj, \pf + k\f). 
We use the same pole dominance idea that we used for the photon propagation terms. In 
this case, the four-vector W appears instead of the four vector k. Then we handle the term 
2pf ■ k' = 2pf ■ W — 2m 2 by using 

' ' + „ J => Pf<W m -W) - 2 - 2 . (59) 



Pf ■ W — m 2 pf ■ W on — m 2 Pf ■ W on — rn 2 Pf • W — m 2 

This allows us to derive the analogous differential equation for K 2 : 



K ^ ' > = 2 P ^£-^ - 2 Pr wl- m ^ ■(«■--(? + *" + M I KS, r). (60) 
where 



K (h r)= f d ' ik ' p -i$-k>)-r 1 F (\ k ~ k '\) / fil N 

K °^ r) ~ 7(2703 6 2p r k' + k' 2 UJ 2_y, + ie - ^ 

Eq. (|60D is equivalent to the full expression for K2 and also shows how we can make a first 
approximation for K 2 by substituting K on {k, r)/(2pj ■ W on — m 2 ) for K 2 on the right hand 
side. The technique is the same as in previous sections. We immediately see that the second 
term vanishes in the separated atom approximation. 

We may evaluate K on in the separated atom approximation, because this is essentially 
the same integral as J on . The separated atom approximation worked except when qi was 
small. Here the quantity \k — \Pf + k\ f\ plays the same role as q\. The result is 

K on (k, r ) = e-^e-^i-^j e W 2E f^ 2+ Pf r F (\k + p f - 'E f u + u 2 + p 2 f\) . (62) 



Since Ef = ^Jm? +pj is very large, F is evaluated with a large argument, and this kills the 
K 2 term. 

The size of the quantity K 2 is controlled by K on and by the denominator 

D = 2p r W on -m 2 = E f {E f +a)- (p) • f) \p f + k\ (63) 

Numerical evaluation shows that K 2 K\. If pf = 0, K on is small because the atomic form 
factor is evaluated at a large argument. When \pf\ takes on a typical experimental value, 
the energy denominator is very large. The net result is that K 2 is ignorable. 

27 



However when \pj\ is very much larger than the electron mass, D approaches and K 2 
can become large. The bremsstrahlung from a collection of atoms would then have a large 
contribution from the electron propagation term. This small value of D is a necessary con- 
dition for the occurance of the Landau-Pomeranchuk-Migdal (LPM) effect JHJ in which the 
long time scale of electron propagation allows a coherent effect which reduces the radiation. 
However, we are concerned with the low energy limit in which the electron momentum is 
much less than its mass. So, for us, the electron propagation term is negligible. 



2. Crossed term 



Suppose we keep the full crossed term. This means that we must evaluate the integral 



d 3 k' 



-i(k—k')-f 



F(\k-k'\) 



(2tt) 3 2 Vl ■ k' - k' 2 + ie k' 2 + ie ^ - k' 2 + ie 
The product of propagators can be written 



(64) 



1 



2 Pi ■ k' - k' 2 k> 2 + ie 



1 1 

+ 



1 



2 Pi ■ k' 



(65) 



_k' 2 + ie 2pi-k'-k' 2 _ 

Again we denote the first term as the photon propagation term and the second term as 
the electron propagation term. The last factor l/(2pi ■ k') vanishes only when the term in 
the bracket vanishes, so this zero residue pole makes no contribution to the integral J2{k, f). 
We may then study two separate integrals 



where 



K 3 (k,r 



J 2 (k,f) = K 3 (k,f) + K i (k,f 



d 3 k' e -^-k')-f l F(\k-k'\) 



(66) 



(2vr) 3 2 Pl ■ k< k' 2 + ie a;?- k' 2 + 



(67) 



and 



K 4 (k, r 



d 3 k' 



-i(k—k')-r 



F(\k-k'\) 



(2tt) 3 2 Pi ■ k' - k' 2 + te 2 Pi ■ k' u 2 - k' 2 + ie 
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(68) 



It is clear that we can handle K 3 using the same techniques that we used for K\. We 
derive the following differential equation for K 3 : 



K 3 (k,r 



Jan{k-i T 



-2pi 



K ■ k 



K 3 {k,r 



(69) 



2pi ■ K 2pi ■ K 

Once again the first term is the separated atom approximation for this particular term. 
This dominates K 3 . For general purposes, it is useful to note that Eq. ( |BTS| ) has the formal 
solution 

1 



K 3 (k,r) 



-J n(k, r ). 



(70) 



2pi ■ k + 2pi ■ (k — k + iV r 
Again, this formal solution gives us a controlled way to study the corrections to the separated 
atom approximation. We use 



with 



(2 Pi ■ K)K 3 (k, r) = J on (k, r) + 6JZ{k, r) 



^J C on{^^) = -2p%- K + iV r -k J on {k,r 



(71) 



(72) 



We find the cross term correction effects due to SJ^ are very small, as illustrated in Fig. 9. 

The term K4 of Eq. (p8[) represents new physics occuring in this two-atom process. We 
use the same techniques we used for K2. There is a pole in the integrand representing the 
long distance propagation of the electron of four-momentum X on = (Ei — u>, \pi — k\f). This 
allows us to derive the differential equation 



K 4 {k,f 
where 
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(2tt)3 2pi ■ k' - k' uj 2 - k' 2 + ie 

Note the appearance of the +m 2 term in the denominator of the first term of K 4 . This 
renders the energy denominator very large, it never vanishes even for infinitely large Ei. 
Careful numerical evaluation leads to negligibly small results for K4. These results are too 
small to be plotted. 
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D. Summary of Bremsstrahlung Cross Section with Thompson Scattering 



The basic expression for the bremsstrahlung holography cross section where the object 
atom scattering can be described by the Thompson amplitude is given by Eqs. fl2~T|) and 
(^2|). The full on- and off-shell integral I(k, f) of Eq. (^) is very well approximated by the 
on-shell integral I on (k,f) of Eq. (f2T)|). 

Furthermore, k^k^ = so that we may perform the sum over the polarization vectors e 
(see Ref. E3, Eq. (7.61)) with the results 
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The net result is that 
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Pf kpi • k p r npi ■ k pfkpf-n pi - kp, L ■ k_ 

The factors in the square brackets account for the peaking of the bremsstrahlung radiation 
intensity which occurs in the direction of the initial electron velocity. This feature represents 
the influence of the vector nature of the photon and is significantly different than the simple 
classical result from scalar electrodynamics. 

The polarization dependent cross section is obtained simply from Eq. (|22|) . by replacing 
the term I(k, f) with I on of Eq. (|26[) . That is, we find 
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VI. BREMSSTRAHLUNG X-RAY HOLOGRAPHY (BXH) INCLUDING 

RESONANT SCATTERING 



We now return to the case where the scattering of the bremsstrahlung photon by the 
object atoms includes the resonance correction terms given by Eqs. (|15D, (|TT|), and (|18D. 
We simply use those equations in the expression for the amplitude given by Eq. (|H]). The 
total Thompson plus resonant hologram amplitude A4t+r is the sum of the Born approxi- 
mation reference term 71 and the Thompson plus resonant object amplitude O. Calculating 



this amplitude and squaring leads to the resonance correction |57W_r| 2 to the non-resonant 



Thompson squared matrix element \A4\ 2 with 



\q | z J {ztt)* 
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\q'\ 2 co 2 -k' 2 + ie 



(79) 

where S(uj) is given by Eq. (|T8"D. If the resonant scattering and the Thompson scattering 
terms are of comparable strength, the full square of the resonant and non-resonant hologram 



amplitude is the sum of the terms given by Eqs. (|77|) and ([79]): 



\M T+R \ 2 = \M\ 2 +\5M R \ 2 . (80) 

The integral over d 3 k' is the same as that evaluated in the previous sections. 
The essential result is that the integral can be evaluated by removing the expression 
V x {k') (1 - F(\q'\)) /\q '| 2 and evaluating it for k! = k. Recall that k = (cu, cur). In that 
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(81) 



As explained above, the diagrams with propagating electrons are also potentially 
troublesome — at high energies these effects produce the interesting LPM effects. Here the 
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off-shell electron effects are governed by known atomic wave functions which go into com- 
puting the function S(u). 

The sum over the polarization vectors of the photon leads to the expression 
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VII. X-RAY FLUORESCENCE HOLOGRAPHY (XFH) 

We now briefly consider the physics of x-ray fluorescence holography (XFH). In this case, 
an atom in the sample is excited from its ground state into an excited state by an incoming 
photon or an incoming electron. After the ionization, the excited atom decays and we must 
consider the interference effects for the outgoing fluorescence photon via the direct path to 
the detector and via the single-scattering paths to the detector. 

The Feynman diagrams for XFH are shown in Fig. 10. Suppose the incident photon 
(Fig. 10a) — or the incident electron (Fig. 10b) — interacts with an atom, knocking an s-shell 
electron into a continuum wave function c. A p-shell electron can spontaneously decay to 
the s-state, emitting a fluorescence photon with the characteristic energies of the atom u = 
E p — E s , and with the reference amplitude 1Z ps f ^ c ^(k, e) for photostimulated fluorescence, 
or with the reference amplitude 1Z es f (c->i)(&j £) for electron stimulated fluorescence. 

First, consider the ionization process. For photoionization, the incoming photon is real or 
on-shell. For electron induced ionization, the incoming electron is on-shell and the ionization 
occurs via the virtual photon exchange, but the potential off-shell photon effects are exactly 
the same as the bremsstrahlung case we have already analyzed in detail. For the relatively 
low energies used in the current XFH experiments, these effects are completely negligible. 

Second, consider the intermediate states. What are the possible off-shell electron effects? 
The electron promoted into the continuum is detectable and is therefore on-shell. The 
virtuality of the continuum electron enters only if there is another final state interaction; 
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such effects are of higher order in a and are neglected here. Thus the only possible off-shell 
electron effects come from the virtual intermediate state electron in the s-state. But this is 
governed by well known atomic wavefunctions. 

For the characteristic radiation used in XFH, the long wavelength approximation holds 
and the radiation is dominated by the electric dipole process. This predominately dipole 
character, combined with the fact that the initial and final electron is on-shell, indicates 
immediately that the separated atom approximation will be extremely close to the exact 
quantum electrodynamic solution. Thus the holographic reference amplitude takes the sim- 
ple form 

Kpsf (c— e)=u) T ps f(k, e) (83) 

for photoionization, and 

TZ esf e) = u T esf (k, e) (84) 

for electron ionization, where the T ps f(k, e) and T es f(k, e) factors account for the remainder 
of the atomic matrix element. 

The holographic object amplitude contribution to the total amplitude occurs because 
the photon is scattered coherently by the object atoms. When the Thompson rescattering 
effects dominate, the on-shell approximations for the FXH holographic interference terms 
are given by the expressions: 

M™ f (k, e) = w T psf (k, e) (l - ± J on (k } f )) (85) 

for photon stimulated XFH, and 

M™ f (k,e) = uj T esf (k,e) (l - ^J w {k,r )) (86) 

for electron stimulated XFH. Here the integral J on (k,r) is given by Eqns. ( p7| ) and (|32|). 

The cross section is obtained by squaring the amplitude and summing over the polariza- 
tion vectors of the photon. Thus the cross section for the intensity of the XFH hologram in 
the on-shell approximation is given generically by 
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^- = u 2 \T(k,e)\ 2 fl-2-ReJ on (k,f)). (87) 
ail \ m / 

And again, when the high electron density regions of the atoms are sufficiently well separated, 
we will recover the usual classical holographic form in the far field limit via Eq. (J3Q[) . 

VIII. MULTIPLE ENERGY X-RAY HOLOGRAPHY (MEXH) 

Finally, we consider very briefly the physics of multiple energy x-ray holography (MEXH). 
In this case, a real photon is sent into the sample from outside and we must consider the 
interference between the direct path to the detector atom and the single scattering paths 
via the object atoms to the detector atom. The Feynman diagrams for this interference are 
shown in Fig. 11. Note that the Feynman diagrams for MEXH are not just the time reversed 
Feynman diagrams for photon induced XFH, and that there are three interfering terms in 
MEXH. 

Within classical electrodynamics, MEXH has been related to XFH by the reciprocity 
theorem, which can be paraphrased roughly as follows: Put the source outside the sample 
and the detector inside the sample, turn on the source, and measure the electric field at 
the detector; if the positions of the source and the detector are interchanged, the electric 
field measured at the detector will be the same. This result comes from the time reversal 
invariance of Maxwell's equations. 

How does this very reasonable classical result emerge from the quantum electrodynamic 
treatment? It clearly is not just simple time reversal invariance, since there are three interfer- 
ing diagrams in MEXH and only two interfering diagrams in XFH. There are three diagrams 
in MEXH because the incoming photons interfere to produce the atomic excitations that 
lead to the fluorescence, and the outgoing fluorescence photons interfere just as they do in 
XFH. However since MEXH must average over many outgoing directions to increase the 



signal level, the interference effects in the outgoing photons will be washed out |32j , and we 
need only consider the first two diagrams in Fig. 11. Then our question becomes how are 
these two diagrams related to the analogous diagrams for XFH shown in Fig. 10. They still 
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are not just the simple time reversed diagrams: in MEXH the Thompson process (or, in 
general, the Compton process) occurs in the incoming state of the photon that will produce 
the photoionization, whereas in XFH the Thompson process is in the outgoing photon state 
of energy (E p — E s ) that will be detected. Since the incoming photon energy used in MEXH 
is not equal to (E p — E s ), the two processes are not related by time reversal invariance. 

We discuss this further by displaying the relevant equation. In MEXH, the holographic 
object amplitude contribution to the total amplitude occurs because the initial photon in the 
incoming beam with momentum kb and polarization e& is scattered elastically by the object 
atoms prior to absorption by the detector atom. If the total matrix element to produce the 
outgoing angle averaged MEXH fluorescence decay is denoted TZ me ^^(kb, 4) an d if the 
Thompson scattering effects dominate, then the on-shell approximation for the holographic 
interference term is given by the expression 

MZ(h,e b ) = (l- -J n(h,f)) u T me (h,e b ), (88) 

where uT me (kb,€b) represents the photon absorption process, and once again the integral 
Jon{kb,r) is given by Eqns. (|27| ) and (|32|). The amplitude for MEXH given by Eq. ( |88| ) 
and the amplitude for FXH given by Eqs. fl8"5] ) and ( |5ED are not complex conjugates of one 
another. However, the physics of these two amplitudes is closely related. 

IX. SUMMARY AND DISCUSSION 

We have shown that, if the energy of the bremsstrahlung photons is measured, the 
bremsstrahlung radiation produced inside a crystal will produce a holographic interference 
pattern in the far field outside the crystal. To use this new form of generalized holography 
f33| , we must know the reference and object amplitudes. These amplitudes were calcu- 
lated using quantum electrodynamics, and compared with the corresponding predictions of 
classical scalar electrodynamics. The essential results for bremsstrahlung holography are 
displayed in Eqs. (|77|), (|7§), and (0). 
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The total amplitude M. is the sum of the reference amplitude 1Z and the object amplitude 
O. To obtain very accurate results, the full expression for the quantity J on (fc, r ) given by 



Eq. (^) must be used. Its simpler asymptotic form given by Eq. ( |30| ) is not accurate for 
the case in which the separation r between the source atom and the object atom is parallel 
to the direction k of the detected photon, as shown in Fig. 5. 

The key feature in obtaining Eqs. (|77f) and fl82f) is that the photons that propagate 
from the source atom to the object atoms are essentially on-shell — the square of their four 
momenta is very close to zero. Sect. V is devoted to the detailed arguments for our on-shell 
separated atom approximation. We show explicitly that all of the known short-ranged off- 
shell virtual effects are negligibly small for the low (40 — 60 keV) electron energies used in the 
current experiments. The underlying reason for this is that the atoms in solids are too far 
apart for the off-shell photons or electrons, produced via bremsstrahlung or via fluorescence, 
to propagate from one atom to another. 

It is interesting to compare the present case in which photons propagate between atoms 
with two examples from nuclear physics: (1) hadronic scattering from nuclei, and (2) pion 
production in nucleon-nucleon or nucleon-nucleus collisions. Beg's theorem applies to 
hadron-nuclear scattering and states that, if the target nucleons are separated by distances 
greater than the range of the hadron-nucleon interaction, then the hadron-nucleus scattering 
amplitude can be expressed in terms of on-shell hadron nucleon scattering amplitudes. This 
is called the separated scatterer approximation [|l^,|l^,|19|,^(| . In this language, our results 
can be stated as the confirmation that the analogous separated atom approximation is valid. 

In our first example, we want to consider the scattering of hadrons from the nucleons 
inside the nucleus. In our condensed matter physics example, we were able to consider 
the scattering of photons from essentially stationary atoms and slowly moving electrons. 
However, in nuclei, the nucleons move extremely rapidly so that the separations between 
the nucleons fluctuate and, in addition, the nucleons can overlap. The average separation 
distance between nucleons is about 1.8 fm, which is about twice as big as the typical range of 
hadron-nucleon interactions (circa 1 fm). However, in the case of very high energy hadronic 
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beams, we can be reasonably sure that the nucleons will not move during the passage of the 
hadron through the target, and that the use of on-shell hadron nucleon amplitudes based 
on the average separation distance is valid 

In our second example, the pion (it) production reactions pp — > dn or p + A — > (v4+l)7r + , 
are the strong interaction analogs of bremsstrahlung holography. In this case, Fig. 2 also 
applies, but the wiggly lines represent pions and the solid vertical lines represent nucleons 
or nuclei. These processes involve high momentum transfer (for pions produced with low 
energy) and small inter-nucleon separations. Therefore off-shell pions can and do propagate 
between nucleons. As a result, no cross sections can be computed to better than about a 



factor of two p5 |. In contrast, in the present paper, we can calculate cross sections reliably 
at the 10 -4 level or better because the atoms and the electrons in the atoms move relatively 
slowly. 

The separability argument works extremely well for the three x-ray holographies (FXH, 
MEXH and BXH) currently under experimental development. In a sense, it is quantum 
electrodynamics that requires the separated atom approximation to work so well since it 
supplies the forces responsible for the relatively slow motion of the atoms and of the electrons 
in the atoms, and it also supplies the interactions between the atoms and the incoming and 
outgoing photons and electrons. 

Thus it is quantum electrodynamics that gives us our main results for bremsstrahlung 
holography summarized in Eqs. (|2lD through (|25|) , Eq. (|77|), and Eq. (|S2"D. Quantum elec- 
trodynamics shows that bremsstrahlung holography should work — the remaining problems 
are experimental. 
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Figure Captions 

1. The classical scalar wave description of internal source holography. The source atom 
produces the spherical reference wave R which propagates directly to the detector, and which 
is scattered by the object atom to produce the spherical object wave O. The interference 
between R and O at the far field detector produces the internal source hologram. 

2. The uncrossed (a) and crossed (b) Feynman diagrams for the reference amplitude TZ 
in bremsstrahlung x-ray holography (BXH). The x represents the target atom that produces 
the bremsstrahlung photon. In the uncrossed (crossed) diagram, the outgoing photon is 
created after (before) the virtual photon is destroyed. The solid line represents the incident 
electron and the wiggly lines the photons; this standard convention is used in all the Feynman 
diagrams in this paper. 

3. The four Feynman diagrams for the BXH object amplitude O. The object amplitude 
has four terms due to the crossed and uncrossed source terms, and the crossed and uncrossed 
Compton scattering terms. The sum of the two 1Z diagrams in Fig. 2 interfere with the sum 
of the four O diagrams in this fi gure to produce the bremsstrahlung hologram. 

4. The electron density p(s) for isolated copper atoms (dashed line) and for copper atoms 
in crystalline copper (solid line) calculated using the FEFF computer code f30fl . Here s is 
the distance from the center of the atom. 

5. Comparison of the real part of the on-shell separated atom approximation J on (solid 
line) given by Eq. ( p7|) with the real part of the classical spherical wave holography function 
(dashed line) given by Eq. (|30|). The results in Figs. 5-9 are shown for crystalline copper 
with representative experimental kinematics: the incident electron energy is 60 keV and the 
outgoing photon energy is 20 keV. Here cos(6 l ) = k ■ f. 

6. The small effects of photon virtuality on the screening correction given by the ratio of 
the second to the first term in Eq. Q4"0]). The ratio Re 5I S /Re [(1 — F(q 2 )) J on ] is plotted to 
illustrate the size of these corrections for two typical experimental values of the momentum 
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transfer, namely 12.2 keV and 93.5 keV. Here the momentum transfer A = p^ — pf and the 
angle 6' is specified by cos{6') = A • f. 

7. The small effects of virtual photon propagation on the Coulomb correction. The 
real part of the on-shell Coulomb correction SJ£° ul /(pf — Pi + k) 2 (dashed line) given by 
Eq. (]4~6D is compared with the real part of the full on-shell separated atom approximation 
Jon/(Pf—Pi + K) 2 (solid line) given by Eq. (p7|) . Here A = 12 keV, jvf = 0.5, and pfPi—0.5. 

8. The small effects of virtual electron propagation in the uncrossed graph in Fig. 2. 
The real part of the uncrossed correction 5J^ cr (dashed line) given by Eq. fl57|) is compared 
with the real part of the on-shell separated atom approximation J on (solid line) given by 
Eq. (pTp. Here k ■ pf — 0.5 and Pi ■ r — 0.5. 

9. The small effects of virtual electron propagation in the crossed graph in Fig. 2. The 
real part of the crossed correction SJ^ (dashed line) given by Eq. ( |72|) is compared with the 
real part of the on-shell separated atom approximation J on (solid line) given by Eq. (|27|). 
Here pi ■ f = 0.5. 

10. The Feynman diagrams for x-ray fluorescence holography (XFH). The s-state core 
hole can be made by photoionization (a), or by electron induced ionization (b). The black 
dot represents the photon-object atom scattering amplitude, c the continuum electron, p 
the p-state electron, and s the s-state electron. Note that the black dot is shorthand for two 
diagrams — namely, the crossed and uncrossed Compton diagrams shown in Fig. 3. 

11. The Feynman diagrams for multiple energy x-ray holography (MEXH). The notation 
is the same as Fig. 10. 
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